Alternative fidelity measure between quantum states 
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We propose an alternative fidelity measure (namely, a measure of the degree of similarity) between 
quantum states and benchmark it against a number of properties of the standard Uhlmann-Jozsa 
fidelity. This measure is a simple function of the linear entropy and the Hilbert-Schmidt inner 
product between the given states and is thus, in comparison, not as computationally demanding. It 
also features several remarkable properties such as being jointly concave and satisfying all of Jozsa's 
axioms. The trade-off, however, is that it is supermultiplicative and does not behave monotonically 
under quantum operations. In addition, new metrics for the space of density matrices are identified 
and the joint concavity of the Uhlmann-Jozsa fidelity for qubit states is established. 
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I. INTRODUCTION 



The understanding of the set of density matrices as 
a Riemannian manifold implies that a notion of dis- 
tance can be assigned to any pair of quantum states. In 
quantum information science, for instance, distance mea- 
sures between quantum states have proved to be useful 
resources in approaching a number of fundamental pro- 
blems such as quantifying entanglement ^ ||, the de- 
sign of optimized strategies for quantum contro l H , ^J, 
and quantum error correction || [?], ||, [)], [h], [0], |T!jf . In 
addition, the concept of distinguishability between quan- 
tum states [^3| can be made mathematically rigorous and 
physically insightful thanks to the close relationship be- 
tween certain metrics for the space of density matrices 
and the error probability arising from various versions of 
the quantum hypothesis-testing problem Q. Distance 
measures are also regularly used in the laboratory to ver- 
ify the quality of the produced quantum states. 

A widely used distance measure in the current litera- 
ture (or, more precisely, a fidelity measure — that is, a 
measure of the degree of similarity — between two gen- 
eral density matrices), is the so-called Uhlmann-Jozsa fi- 
delity, 3\ Historically, this measure had its origins in 
the 1970s through a set of works by Uhlmann and Al- 
berti Jl5|, [I?], |l8) , who studied the problem of gener- 
alizing the quantum mechanical transition probability to 
the broader context of ^-algebras. The use of the term 
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fidelity to designate Uhlmann's transition probability for- 
mula is much more recent and initiated in the works of 
Schumacher [[l9| and Jozsa [^Oj. Indeed, in an attempt 
to quantify the degree of similarity between a certain 
mixed state p and a pure state Schumacher dubbed 
the transition probability (ip\p\ip) the fidelity between 
the two states. In parallel, Jozsa recognized Uhlmann's 
transition probability formula as a sensible extension of 
Schumacher's fidelity, where now the measure of similar- 
ity is related to a pair of mixed states p and cr. Ever 
since, Uhlmann's transition probability formula has been 
widely accepted as the generalization of Schumacher's fi- 
delity. 

The prevalence of this measure as one of the most used 
notions of distance in quantum information is not acci- 
dental, but largely supported on a number of required 
and desired properties for the role. For example, 3" sa- 
tisfies all of Jozsa's axioms, that is, besides recovering 
Schumacher's fidelity in the case where one of the states is 
pure, the following three additional properties also hold: 
First, "5 equals unity if and only if it is applied to two 
identical states; in other cases it lies between and 1. 
Second, it is symmetric, i.e., the fidelity between p and 
cr is the same as that between a and p. Third, it is 
invariant under any unitary transformation on the state 
space. Nevertheless, 1 is not the unique measure satis- 
fying these properties. A prominent alternative which 
also complies with Jozsa's axioms and shares many other 
properties of 5F, is given by the nonlogarithmic variety of 
the quantum Clrcrnoff bound, Q, recently determined in 
Ref. |2l| . In analogy with its classical counterpart p2| , 
the quantum Clrcrnoff bound determines — in the limit 
of asymptotically many copies — the minimum error pro- 
bability incurred in discriminating between two quantum 
states |2l], |2^| . 

Despite fulfilling the properties listed above, both 
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and Q are, in general, unsatisfying measures from a prac- 
tical computational viewpoint. Although 3 r can be ex- 
pressed in a closed form in terms of p and <x, it involves 
successive computation of the square roots of Hermi- 
tian matrices, which often compromises its use in ana- 
lytical computations and numerical experiments, espe- 
cially when the fidelity measure must be computed many 
times. Even more serious is the case of Q, which to date 
has only been defined variationally as the result of an 
optimization problem p3| . The question that naturally 
arises is whether an easy-to-compute generalization of 
Schumacher's fidelity can be obtained. In this paper, we 
provide a positive answer to this question and a thorough 
analysis of our proposed alternative fidelity measure 3^ jy . 

Recently, we became aware of the very recent work of 
Miszczak et al. pij in which 3> was introduced as an up- 
per bound to the Uhlmann-Jozsa fidelity. In many ways 
our analysis of 3" N is complimentary to that provided in 
Ref. [p4[ ; results in common are noted in the correspon- 
ding sections of our paper. 

Our paper is structured as follows. In order to provide 
a concrete ground for our proposal of 3" N as an alterna- 
tive fidelity measure, we first reexamine, in Sec. H a set 
of basic properties of the Uhlmann-Jozsa fidelity. In Sec. 
[II we formally introduce 3^ and analyze it in the spirit 
of the properties reviewed in Sec. [n|. The computational 
efficiency of 3^ is contrasted with a number of previ- 



IV. We 



ously known distance or fidelity measures in Sec. 
summarize our main results and discuss some possible 
avenues for future research in Sec. [v| 



II. UHLMANN-JOZSA FIDELITY 

In this section, we will briefly survey some physically 
appealing features inherent to the Uhlmann-Jozsa fidelity 
3". In Sec. Ill, these features will be used as a reference 



for characterizing the proposed fidelity measure. 



A. Preliminaries 

The Uhlmann-Jozsa fidelity 3 r was originally intro- 
duced as a transition probability between two generic 
quantum states p and cr ]15|: 



?(p,a) :=max \{^\f)\ 2 = 



Tr JJptTy/p 



(1) 



Here, \ip) and \tp) are restricted to be purifications of 
p and cr, while the second equality indicates that the 
maximization procedure can be explicitly evaluated. At 
this stage, it is worth noting that it is not uncommon to 
find Vj being referred to, instead, as the fidelity (see, 
for example, Ref. [p2|). 

In Ref. |2(|, Jozsa conjectured that Eq. ([![) was the 
unique expression that satisfies a number of natural pro- 
perties expected for any generalized notion of fidelity p5| . 
Throughout, we shall refer to these as Jozsa's axioms: 



1. Normalization — i.e., ^(p, cr) € [0, 1] with the up- 
per bound attained iff p = cr (the identity of indis- 
cernible property). 

2. Symmetry under swapping of the two states — i.e., 
?(p,cr)=2(cr,p). 

3. Invariance under any unitary transformation U of 
the state space — i.e., 3 r (XJpXJ' f , UerUt) = { J{p, cr). 

4. Consistency with Schumacher's fidelity when one of 
the states is pure — i.e., 



?(p,|VXV>|) = <V'IpI 

for arbitrary p and \ip). 



(2) 



The proof that 3 r satisfies all of Jozsa's axioms follows 
easily from the variational definition of Eq. (Q) (see, e.g., 
Ref. ^] for technical details). The remainder of this 
section discusses a number of less immediate properties 
of 3\ 



B. Concavity properties 

The concavity property of quantities like entropy, mu- 
tual information, and fidelity measure are often of the- 
oretical interest in the quantum information community 
p2[ . In this regard, it is worth noting that a useful fea- 
ture of 3" is its separate concavity in each of its arguments; 
i.e., for pi,p2 > 0, p\ +P2 = 1, and arbitrary density ma- 
trices p 1 , p 2 , <Tx, and <x 2 , we have 

?(PlPl +V2pi,cr 1 ) >p 1 J{p 1 ,cr 1 )+p 2 J(p 2 ,cr 1 ). (3) 

By symmetry, concavity in the second argument follows 
from Eq. (^) . Separate concavity can be proven |2(]] 
using the variational definition of 3 r from Eq. (|l|). 

While it is known that \fj is jointly concave |ll| po| , 
i.e., 



+ P2P 2 ,Pl< T l +P2Cr 2 ) 
> p 1 Vj{p 1 ,cr 1 ) +p 2 V / J(p 2 7' T 2), 



(4) 



it is also known that the Uhlmann-Jozsa fidelity 3" does 
not, in general, share the same enhanced concavity prop- 
erty |27[|. 



C. Multiplicativity under tensor products 

Another neat mathematical property of 3 r (p, cr) is that 
it is multiplicative under tensor products: for any density 
matrices p 1 , p 2 , cr 1 and cr 2 : 



?(Pl ® Pi, crx ® cr 2 ) = 3 r {p 1 ,cr 1 )3 : (p 2 ,cr 2 ) 



(5) 



This identity follows easily from the following facts: for 
any Hcrmitian matrices A and B, (i) Tr(A <E> B) = 
Tr(A) Tr(B) and (ii) VA <E> B = \/A ® y/B. 
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An immediate consequence of this result is that for two 
physical systems, described by p and cr, a measure of 
their degree of similarity given by 5" remains unchanged 
even after appending each of them with an uncorrelated 
ancillary state r — i.e., 5'(p ® r, cr (g) r) = 5"(p, cr). 



D. Monotonicity under quantum operations 

Given that 1(p, cr) serves as a kind of measure for the 
degree of similarity between two quantum states p and 
cr, one might expect that a general quantum operation 
£ will make them less distinguishable and, hence, more 
similar according to |22| : 

y(£(p),£(<r))>?(p,o-)- (6) 

Indeed, it is now well known that Eq. (Q) holds true 
for an arbitrary quantum operation described by a com- 
pletely positive trace-preserving (CPTP) map £ : p 
£(p). Inequality (^) qualifies J as a monotonically in- 
creasing measure under CPTP maps and can be consi- 
dered the quantum analog of the classical information- 
processing inequality — which expresses that the amount 
of information should not increase via any information 
processing. 

On a related note, it is worth noting that any mea- 
sure M which is (i) unitarily invariant, (ii) jointly con- 
cave (convex) , and (iii) invariant under the addition of an 
ancillary system is also monotonically increasing (decrea- 
sing) under CPTP maps [p8| . Clearly, since satisfies 
all the above-mentioned conditions, Eq. (||) also follows 
by simply squaring the corresponding monotonicity ine- 
quality for 



F. Trace distance bounds 

An important distance measure in quantum informa- 
tion is the metric induced by the trace norm |j ■ || tr , which 
is commonly referred to as the trace distance p2| : 

D(p,ir) = |||p-tr|| tr . (10) 

The trace distance is an exceedingly successful distance 
measure: it is a metric (as is any distance induced by 
norms), unitarily invariant 63, jointly convex p^[ , de- 
creases under CPTP maps [p7f , and in the qubit case, 
is proportional to the Euclidean distance between the 
Bloch vectors in the Bloch ball. The trace distance is 
also closely related to the minimal probability of error 
on attempts to distinguish between a single copy of two 
nonorthogonal quantum states ]3l| . For all of these rea- 
sons, one is generally interested to determine how other 
distance measures relate with the trace distance. 

The following functions of the Uhlmann-Jozsa fidelity 
were shown in Ref. Q to provide tight bounds for D 
& 

1 - v^pT^) < T>(p, cr) < y/1 - ?(p, cr). (11) 

In fact, the stronger lower bound 1 — 3 < D holds if p and 
cr have support on a common two-dimensional Hilbcrt 
space jll| (e.g., any pair of qubit states) or if at least one 
of the states is pure p3j. 

From these inequalities, one can conclude a type of 
qualitative equivalence between the Uhlmann-Jozsa fi- 
delity "J and the trace distance D: whenever J is small, 
D is large and whenever J is large, CD is small. 

III. ALTERNATIVE FIDELITY MEASURE 
A. Preliminaries 



E. Related metrics 



The Uhlmann-Jozsa fidelity by itself is not a metric 
(for a quick review of metrics, see Appendix |a|) . However, 
one may well expect that a metric, which is a measure 
of distance, can be built up from a measure of similarity 
such as 3\ Indeed, the functionals 



A[3{p, cr)] := arccos ( a/J(p, cr) 



B[?(p,er)\ :=\j2-2y/3{p,cr), 
C[3{p,cr)] := y/1 -5(p,,r) 



(7) 

(8) 
(9) 



exhibit such metric properties (see Refs. [£2[ |3l], |32|, |3^, 
|3~fl . |il5|] and also Appendix [B 3| for more details). In par- 
ticular, these functionals arc now commonly known in 
the literature, respectively, as the Bures angl e [p2| , the 
Bures distance 132, p33|] , and the sine distance |35|. 



We shall now turn attention to our proposed measure 
of the degree of similarity between two quantum states p 
and cr — namely, 

7 s {p, cr) = Tr(pcr) + y/l - Tr(p 2 ) y/l - Tr(cr 2 ). (12) 

This is simply a sum of the Hilbcrt-Schmidt inner pro- 
duct between p and cr and the geometric mean between 
their linear entropies. It is worth noting that the same 
quantity — by the name superfidelity — has been inde- 
pendently introduced in Ref. p4| as an upper bound for 
J. 

Remarkably, when applied to qubit states, S N is pre- 
cisely the same as 3\ This observation follows easily from 
the fact that for density matrices of dimension d = 2, it 
is valid to write 

3 N (p, cr)| rf=2 = Tr(p<x) + 2v/dct (p)Vdct (cr), (13) 

which is just an alternative expression of !F for qubit 
states fl3l El. 
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When d > 2, however, 3> no longer recovers 5F, but 
can be seen as a simplified version of the fidelity measure 
proposed by Chen and collaborators |43| , which reads as 



l-r 1 



(14) 



where r = l/{d— 1) and d is the dimension of the state 
space of p and er. Moreover, it is straightforward to verify 
that while 3> reduces to the Schumacher's fidelity [the 
right-hand side of Eq. (^)] when one of the states is pure; 
the same cannot be said for £F . 

It is not difficult to see from Eq. ( |l2| ) that 3 N sa- 
tisfies Jozsa's axioms 2, 3, and 4 as enumerated in Sec. 
[I A. The non-negativity of J N required by axiom 1 is 
also immediate from the definition. As a result, 3> is 
an acceptable generalization of Schumacher's fidelity ac- 
cording to Jozsa's axioms if the following proposition is 
true. 

Proposition III.l. 3 N (p,cr) < 1 holds for arbitrary 
density matrices p and u , with saturation if and only if 

P = (T. 

Proof. To begin with, recall that any dxd density matrix 
can be expanded in terms of an orthonormal basis of 

j2_-| 

Hermitian matrices {vk}i -— n such that Tr(viVj) = 5ij 
(see, for example, Refs. JS], (45[). In particular, if we let 
T := (vq, . . . ,Vcp_x), then p and er admit the following 
decomposition: 



p = f ■ T and a = s ■ T, 



(15) 



where f and s are real vectors with d 2 entries (corres- 
ponding to the expansion coefficients which can be de- 
termined using the orthonormality condition). Since p 
and er are density matrices, rand s satisfy < r ■ s < 1 
and r, s < 1, where r = \\r\\ and s = ||s||. 

Using the expansion of Eq. ([Hi]) in Eq. (|l2|), we arrive 
at the following alternative expression of 5F N , 



f N (r, s) = f ■ s + \f\ 
= R-S, 



(16) 
(17) 



where, in the second line, we have defined two unit vec- 
tors in M. d +1 , explicitly, 



R := ( r, \J\ -r 2 ) and S := (s, v 7 ! - s 2 



(18) 



The normalization of R and S implies that 3 r N (p, <r) = 
R ■ S < 1, with saturation if and only if R = S, or 
cquivalcntly p = a. □ 



B. Concavity properties 

As with y/W, the measure ^ N is jointly concave in its 
two arguments; i.e., for pi,p2 > 0, p± + P2 = 1, and 



arbitrary density matrices p 1; p 2 ; an d <T 2, we have 

9w(PlPl +P2P 2 i-Pl cr l +P20"2) > 



pi3>(p 1 ,cr 1 ) +p 2 3 r jv(p 2 > <T 2) 



(19) 



Since ? fails to be jointly concave in general, has 
a stronger concavity property. Remarkably given the 
equivalence between 5" and 3> in the d = 2 case, the 
result of this section implies that !F is jointly concave 
when restricted to qubit states. 

The rest of this section concerns a proof of this con- 
cavity property of 9> . We start by proving the following 
lemma, which provides a useful alternative expression of 
inequality (|l9|). 

Lemma III.l. Define a function F : [0, 1] — > M fry 

F(x) := (r + xu) ■ (s + xv) 



+ VI - ||f + V 1 - P + ^ll 2 - (20) 
Given the density matrices p l7 p 2 , <J\, and cr 27 t/iere 



exisi vectors r, s, u, v € 
inequality 



b el- 



and x £ [0, ll such that the 



F(x)>(l-x)F(0)+xF(l) (21) 
is equivalent to Eq. (|l9|). 

Proof. The proof is by construction. Using the parame- 
trization of Eq. ([l5]) for the density matrices in inequality 
(|l9|), we obtain the following equivalent inequality for the 
vectors rv and s^: 



f N {piri +p 2 r 2 ,pisi +P2S2) > 
Pi/jv(^i,si) +p 2 /jv(r 2 , s 2 ), 



(22) 



where the function f N was defined in Eq. (|l6|). 

A straightforward computation shows that inequality 
( pl| ) is identical to inequality (|22|) when we identify x = 
P2, 1 — x = px, and set 



r = T\ , u = Ti — r\ , 

S = S 1 , ?7 = S2 — S; 1 . 



(23) 
□ 



If f (a:) has negative concavity in a; G [0, 1], then ine- 
quality ( |2l| ) is automatically satisfied as it establishes 
that the straight line connecting the points (0, F(0)) and 
(1,F(1)) lies below the curve {(x,F(x))\x £ [0,1]}. As 
a result, the joint concavity of 3 N is proved with the 
following proposition. 

Proposition III. 2. For x G [0, 1], and r, s, u, v G M. d " 
specified in Eq. (jjg|), i/ie function F(x) [cf. Eq. 
satisfies 



d 2 F{x) 
dx 2 



< 



(24) 



and hence 3 N is jointly concave. 



The proof of this proposition is given in Appendix B 1 
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C. Multiplicativity under tensor product 

In contrast with 5F, the new measure J N is not multi- 
plicative under tensor products. In fact, it is generally 
not even invariant under the addition of an uncorrelated 
ancilla prepared in the state r. In this case, 3> between 
the resulting states reads as 

3 N (p®T,cr®T) = Tr(pa) Tr(r 2 ) 



for the resulting pair of states? An affirmative answer 
would follow from a proof of the inequality 



+ y/l- Tr(p 2 ) Tr(T 2 )v/l - Tr(<x 2 ) Tr(r 2 ), 

where the left-hand side equals 3 N (p, er) if and only if 
Tr(r 2 ) = 1 or, in other words, if and only if r is a pure 
state. More generally, it can be shown that 3> is super- 
multiplicative, i.e., 

3«(Pi <8> P 2 , (T i ® er 2 ) > 3 : N (p 1 ,(Ti)J N (p 2l (T2). (25) 

A proof of this property is given in Appendix |B 2[ a si- 
milar proof was independently obtained in Rcf. 



D. Monotonicity under quantum operations 

That 3 r N is only supermultiplicativc may be a first sign 
that it may not behave monotonically under CPTP maps. 
In fact, as we shall see below, Ozawa's counterexample 
p6| to the claimed monotonicity of the Hilbert-Schmidt 
distance [f|7) can also be used to show that J N does not 
behave monotonically under CPTP maps. 

Let p and <x be two two-qubit density matrices, written 
in the product basis as 



and 





(26) 



and consider the (trace-preserving) quantum operations 
of tracing over the first or second qubit. A straightfor- 
ward computation shows that if the first qubit is traced 
over, then 

J^TnO^Tntar)) = 1 > i = J N (p,a), (27) 

which satisfies the desired monotonicity property. Ho- 
wever, if instead the second subsystem is discarded, we 
find 



J N (Tr 2 (p),Tr 2 (ar)) = < i = $ N {p,Z). 



(28) 



Together, Eqs. (27) and (|28| ) show that J N is neither 
monotonically increasing nor decreasing under general 
CPTP maps. 

A natural question that follows is whether 3 N features 
a weaker form of monotonicity. For example, do arbi- 
trary projective measurements — with the measurement 
outcomes forgotten — give rise to a higher value of 3 N 



(29) 



for any complete set of orthonormal projectors P, and 
for arbitrary density matrices p and er. 



It is a simple exercise to prove Eq. (29) for the particu- 
lar case where either of the commutation rules [Pi , p] = 
or [P^ er] = is observed for all values of i. Whether the 
same conclusion can be drawn for the more general, non- 
commutative cases remains to be seen. In this regard, 
we note that a preliminary numerical search favors the 
validity of Eq. (§|). 



E. Related metrics 

In parallel to th e metrics A[5F], B[3], and Cf?] intro- 
duced in Sec. HE, we define 



A[3 f N (p,a)] := arccos ( yJ$ N {p, ctH , (30) 

B[5 N {p,a)\ := A /2-2V^(pT^y, (31) 
C[9 N (p, cr)}:= y/l - 5„(p,<r), (32) 

and prove that while Cp^] preserves the metric pro- 
perties, both Apjy] and do not always obey the 
triangle inequality 

X[J n ( P ,ct)]<X[J n {p,t)]+X[^ n (t,(t)], (33) 

where X here refers to either A, B, or C . For example, 
consider the qutrit density matrices, p = I3/3, 



0.90 0.04 0.03 
0.04 0.05 0.02 I . (34) 
0.03 0.02 0.05 



TABLE I: A numerical test of the triangle inequality for 
A[? N ], B[? N ], and C[J N ]. 








°\ 1 


I 





I and t = 1 










X 


X[fs(p,*)] 


X[3 n ( P ,t)] + X[J n (t,<t)] 


A 


0.9553 


0.9241 


B 


0.9194 


0.9137 


C 


0.8165 


0.8828 



Numerical computation of the quantities appearing in 
the triangle inequality gives rise to Table H. Note that for 
X = A, B, the first column dominates the second; i.e., 
the triangle inequality is violated and therefore neither 
AfJjy] nor ^[Jjv] is a metric. For X = C, no violation is 
observed for the above density matrices. Next, we prove 
that this is the case for any three density matrices p, a 
and t; thus, C[3>f] is a metric. 
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Proposition III. 3. The quantity C[3 r N (p, cr)] is a me- 
tric for the space of density matrices. 

To prove this proposition, we will make use of the 
following theorem due to Schoenberg (see also |4£], 
Chap. 3, Proposition 3.2]). We state here an abbreviated 
form of the theorem sufficient for our present purposes. 

Theorem III.l (Schoenberg). Let X be a nonempty set 
and K : X x X — > M a function such that K(x, y) = 
K(y, x) and K(x, y) > with saturation iff x = y, for all 
x, y € X. If the implication 



= 



K{xi, Xj) CiCj < 



(35) 



holds for all n > 2, {xi, . 
K, then \f~K is a metric. 



,x n } C X, and {ci, . . . , c„} C 



We make a small digression at this point to remark 
that, in spite of its successful application on the grounds 



51 



of classical probability distance measures |5C 
Schoenberg's theorem has received almost no attention 
by the quantum information community. In this paper, 
besides proving the metric properties of C^n], we will 
also make use of Schoenberg's theorem to provide inde- 
pendent proofs of the metr ic p roperties of B[3 r (p, cr)] and 
C[1(p, cr)] (see Appendix ||J). 



Proof of Proposition III. j| . Clearly, from the definition 
of C 2 [9 r JV (p, cr)], it is easy to see that it inherits from 
3 N (p, cr) the property of being symmetric in its two ar- 
guments and that C 2 [J N (p, c r)] > with saturation iff 
p = cr. So, to apply Theorem [II. 1, we just have to show 
that for any set of density matrices {Pi}™ =1 (n > 2) and 
real numbers {ci}™ =1 such that X)"=i c « = ®> ^ is true 
that 



53 p J] * Cj <0. 

i,j=l 



(36) 



This follows straightforwardly by exploiting the zero-sum 
property of the (real) coefficients Cj and the linearity of 
the trace, 



£ 

: - Tr 



1 - Tr(p iP ,) -Jl- Tr(p?Wl - Tr(p 



53 Cip > 



53c jV /l-Tr(p 2 ) 



< 0, 



which concludes the proof. 



(37) 
□ 



We note that a proof of the metric property of 
V / 2C[9 r JV (p, cr)] — by the name modified Bures distance 
- was independently provided by Ref. j24j. The proof 
provided above is significantly shorter thanks to the 
power of Schoenberg's theorem. 



Trace distance bounds 



In Sec. II F, we have seen that a kind of qualitative 



equivalence between D and 5" can be established through 
the bounds on D given by functions of £F; cf . Eq. (|ll]) . 
Here, we will provide similar bounds on D in terms of 
functions of J N . 

Proposition III. 4. For any two density matrices p and 
cr of dimension d, the trace distance T>(p, cr) satisfies the 
following upper bound: 



D(p,<7) < J-^l-? N (p,(T), 



(38) 



where r := rank(p — cr). Moreover, this upper bound on 
D can be saturated with states of the form 

Udiag[A d ]Ut Udiag[P(A d )]Ut 

P = m , , r . i- and cr = — — , , — — r- — , (39) 

H Tr(diag[A d ]) Tr(diag[A d ]) 1 ; 

where U is an arbitrary unitary matrix of dimension d, 
Ad is an ordered list of d elements taking values in the 
set {Ai,A2} (Ai,A2 > 0, but not simultaneously zero) 
and P(Ad) is the list formed by some permutation of the 
elements in Ad- 

Proof. Note that the product of square roots in the ex- 
pression of 3 Nl Eq. (|T2]), is the geometric mean between 
the linear entropies of p and cr. It then follows from the 
inequality of arithmetic and geometric means that 

1 - Tr(p 2 ) 1 - Tr(<x 2 



> v/l-Tr(p 2 )v/l-Tr(<T2), 

(40) 

which can be reexpressed as the following inequality after 
summation of Tr(pcr) to both sides: 



|p - <t|| hs < y/2[l -9 N (p,tr)]. 



(41) 



Here, ||X|| H s := ^Tr(X^X) is the Hilbert-Schmidt norm 
(also known as Frobenius norm), defined for an arbitrary 
matrix X. The Hilbert-Schmidt norm and the trace norm 

HXlltr 



Tr(VXtX) are related according to 



|X||tr < x/y||x|| 



(42) 



where y := rank(X). Used in Eq. (41), the above inequa- 
lity leads to the desired result 



D(p,er) = f||p- cr|| tr < ^l-y N {p,cr). 



(43) 



To prove that the states in Eq. Q39|) saturate this 
bound, we first note that because those states are isospec- 
tral, their linear entropies are identical and hence inequa- 
lity (ptil ) is saturated. T o prove saturation of inequality 
(p2|), simply use Eq. 



C tr 



Tr 



j) to compute 
r|Ai 



A 2 | 



|p-o-|| HS = VTr[(p-cr)2] 



21 - 



Tr(diag[A d ])' 

Vt|Ai-A 2 | 
Tr(diag[A d ])' 



(44) 
(45) 
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FIG. 1: (Color online) Plot of the trace distance T)(p,a) ver- 
sus 1 — 9 r N (p,cr) for 4 x 10 6 pairs of randomly generated p 
and cr with (a) d = 3 and (b) d — 6. The darker (blue) points 
are generated using pairs of mixed states whereas the lighter 
(green) points are generated using at least one pure state. 
The antidiagonal solid line is the conjectured lower bound, 
whereas the upper bounds given by Eq. (Bq) are represented 
by the dashed curves (cyan) — one for each integer value of 
r 6 [2,4 A gap can be clearly noticed in plot (a) between 
the distribution of states and the absolute upper bound — i.e., 
the right-hand side of inequality (^) with r = d. Such a gap 
occurs whenever d is odd. However, no gap is observed in plot 
(b) between the bulk of randomly generated states and the 
absolute upper bound. In fact, this bound can be saturated by 
density matrices of the form given by Eq. ( [39|) whenever d is 
even. 



from which the identity 
immediate. 



V~t\\ P 



|hs is 

□ 



How good are these upper bounds? With some 
thought, it is not difficult to conclude that the states 
arising from Eq. (39) can only have even r and are thus 
unable to saturate the upper bound of Eq. ( |38| ) for odd 
r. Nonetheless, from our numerical studies, it seems like 



the absolute upper bound — corresponding to the choice 
r = d on the right-hand side of Eq. ( |38| ) — is actually 
unachievable by any states if d is odd. An illustration of 
this peculiarity can be seen in Fig. 1(a), where the up- 
per bound corresponding to r = 3 is well separated from 
the region attainable by physical states. In contrast, for 
every even d, the states given by Eq. (|3^) do trace out 
a tight boundary for the region attainable with physical 
states, as shown in Fig. 1(b) for d = 6. 

On the other hand, it can also be seen from Fig. 1 
that no points occur in the region where T> < 1 — 3>. 
Indeed, intensive numerical studies for d = 3, 4, . . . , 50 
have not revealed a single pair of density matrices which 
contributed to a point in this region. This suggests that 
the following lower bound on D, in terms of 1 N , may well 
be established. 

Conjecture III.l. The trace distance D(p, cr) and the 
measure $ N (p,cr) between two quantum states p and cr 
satisfy 



T>(p,<r) >l-$ N (p,a). 



(46) 



In relation to this, it is also worth noting that the follow- 
ing (weaker) lower bound can readily be established via 
a recent result given in Ref. p4fl : 

Proposition III. 5. The trace distance D(p, cr) and the 
measure r J N {p,cr) between two quantum states p and cr 
satisfy the inequality 



T)(p,a) > 1- y/¥s{p,&). 



(47) 



Proof. This lower bound on D follows immediately from 
the lower bound on D given in inequality (|ll"| ) and the 
inequality 7 < 7 N recently established in Ref. [^4j . □ 

As with the Uhlmann-Jozsa fidelity we can thus 
infer that whenever 7 N is large enough, D is close to zero 
and whenever 5" N is close to zero, D is close to unity. 
However — as should be clear from Fig. 1(b) — the 
converse implication is not necessarily true. 



IV. COMPUTATIONAL EFFICIENCY 

For two general density matrices p and cr, analytical 
evaluation of the Uhlmann-Jozsa fidelity 5"(p, cr) can be a 
formidable task. This is in sharp contrast with 3>(p, cr), 
which involves only products and traces of density matri- 
ces. Even at the numerical level — due to the complica- 
tion involved in evaluating the square root of a Hermitian 
matrix — the computation of J(p, cr) can be rather re- 
source consuming. For a quantitative understanding of 
the computational efficiency, we have performed a nu- 
merical comparison of the time required to calculate the 
fidelity measures S" and 3 N , the nonlogarithmic variety of 
the quantum Chcrnoff bound Q, and the trace distance 
D. We have implemented the computations in both Mat- 
lab and C; we present the Matlab codes for reasons of 
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FIG. 2: (Color online) A semilog plot of the average compu- 
tation time in Matlab and C for the fidelity measures 3" (O)i 
3 N (x), the nonlogarithmic variety of the quantum Chernoff 
bound Q (+), and the trace distance T) (□) as a function of the 
dimension d of the state space. The smaller and larger mark- 
ers correspond to timings from Matlab and C, respectively. 
Computations were performed on a 2.6-GHz Intel Pentium 4 
CPU. 



accessibility and succinctness, while the C codes provide 
accurate timings without the overhead of the Matlab in- 
terpreter. 

The time required to evaluate each function was esti- 
mated by averaging the times for 100 pairs of randomly 
generated <i-dimensional density matrices J54|. Results 
are shown in Fig. 2 as a function of d. The Matlab codes 
arc presented in Appendix |C| we attempted to make 
these codes as efficient as possible within the constraints 
of the Matlab environment. Corresponding C codes were 
implemented as Matlab MEX-filcs for convenience and 
can be found online J56| . Our C implementation di- 
rectly calls the LAPACK and BLAS libraries included 
in the Matlab distribution for eigenvalue decompositions 
and matrix operations. The minimization required in the 
computation of Q was performed using the Brent mini- 
mizer from the GNU Scientific Library J57j . 

The results shown in Fig. 2 are consistent with the 
expected algorithmic complexity: Both J and Q re- 
quire two Hcrmitian diagonalizations, taking an expected 
0(d 3 ) operations each |58|. Computing Q is the slowest 
since it requires both sets of eigenvectors, while C? re- 
quires only eigenvalues from one of the diagonalizations. 
Next fastest is the computation of CD which requires only 
eigenvalues from a single diagonalization. All of 5F, Q, 
and CD appear to require an asymptotic 0(d 3 ) operation 
because of the necessity of diagonalization or some other 
method for computing functions of the input matrices. 



On the other hand, our proposed fidelity measure J N re- 
quires only three Hilbert-Schmidt inner products, with 
asymptotic performance 0(d 2 ). Figure 2 clearly shows 
that the practical numerical evaluation of 3> is dramat- 
ically faster than the evaluation of t J 1 CD, or Q. This 
raises the prospect of using $ N as a numerically efficient 
estimate of distance measures such as and CD - 

particularly for small d where the bounds proven in Sec. 
Ill F arc tighter. As the dimension increases, the com- 
putational advantage of using J N becomes even greater, 
but the quality of the estimate drops. 



V. CONCLUDING REMARKS 

In this paper, we have proposed the quantity 3> as 
an alternative fidelity measure (namely, a measure of the 
degree of similarity) between an arbitrary pair of mixed 
quantum states. This measure, the prevailing Uhlmann- 
Jozsa fidelity J , and the nonlogarithmic variety of the 
quantum Chernoff bound Q |H| are, to the best of our 
knowledge, the only known fidelity measures between 
density matrices that comply with Jozsa's axioms [ p0| . 
That is, 3", Q, and 3 N are the only known measures that 
generalize to pairs of mixed states the concept of fidelity 
introduced by Schumacher between a pure and a mixed 
state pf. 

The simplicity of J N is in sharp contrast with 2f and 
Q since it involves only products of density matrices. 
Numerically, this leads to significant reduction in com- 
putation time for Cf N (p, cr) over 3(p,cr), especially for 
higher-dimensional systems. 

Besides being easier to compute, 3 N has also been 
shown to preserve (and even enhance) a number of the 
useful properties of 2f and Q. For example, we have 
shown that J N is a jointly concave measure, that it can 
be used to place upper and lower bounds on the value of 
the trace distance, and that it gives rise to a new metric 
for the space of density matrices. A remarkable con- 
sequence of the joint concavity of 3 N is that J is also 
jointly concave when restricted to a pair of qubit states 
— an interesting problem which remained unsolved thus 
far §| 0. 

Our measure, nevertheless, is not without its draw- 
backs. To begin with, 5F N — unlike measures such as 
J or Q — does not behave monotonically under CPTP 
maps. In addition, it does not necessarily vanish when 
applied to any pair of mixed states, which are otherwise 
recognized to be completely different according to J, Q, 
or their trace distance CD. In fact, the explicit dependence 
on the linear entropies of p and <r gives rise to the fol- 
lowing undesirable feature: the value of between two 
completely mixed states residing in disjoint subspaces can 
get arbitrarily close to unity as the dimension of the state 
space tends to infinity. 

The undesirable features of J N provide a clue as to 
when Cfjy may not be the preferred measure of similarity 
between two quantum states: We know that J N does not 



9 



measure the similarity between two high-dimensional, 
highly mixed states (i.e., states having non- negligible li- 
near entropy) in the same way that measures like 3, Q, 
or D would. In these cases, the interpretation of 3>, as 
a measure of similarity between quantum states must be 
carried out with extra caution. 

With this in mind, we nevertheless see 3> as an at- 
tractive alternative to 3\ Even when out of its range of 
applicability, it follows from a very recent result of Mis- 
zczak et al. |2^| that 3> provides an upper bound on the 
Uhlmann-Jozsa fidelity 3\ Moreover, it seems promising 
that J N between any two quantum states may be mea- 
sured directly in the laboratory, without resorting to any 
state tomography protocol j24| . 

Let us now briefly mention some possibilities for future 
research that stem from the present work. To begin with, 
it would be interesting to search for a quantitative rela- 
tionship between $ N and Q analogous to that between 
$ N and D established in this paper or that between J N 
and 1 given in Ref . |24j . An estimate of Q based on some 
function of J N would be useful given that a closed form 
for Q is not currently known and that 3> can be com- 
puted relatively easily. In addition, assuming 3~ N as an 
alternative to £F, it seems reasonable to reexamine some 
of the problems where J has proven useful, but with J N 
playing its role. In particular, it would be interesting 
to investigate whether the simplicity associated with 3^ N 
will offer some advantages over 3". 

As a first example, we recall from Refs. ^ ^ that 
a standard measure for the amount of entanglement of 
a state p is given by the shortest distance from p to 
the set of separable density matrices. Given the relative 
simplicity of 3> with respect to 5F, it is not inconceivable 
that a distance measure based on 3 N (such as C[3>]) may 
lead to a more efficient determination of this quantity if 
compared, for example, to C\fi\ or the Bures distance ||. 
Of course, since Cp^] docs not satisfy all the sufficient 
conditions required to give rise to a good entanglement 
measure Q , any serious attempts in this direction should 
be preceded by further investigation of the impact of the 
nonmonotonicity of 3> under CPTP maps. For instance, 
the nonmonotonicity of 3> might also imply that the 
shortest distance from any given state p to the set of 
separable states — as measured by C[3>] — does not 
satisfy the necessary conditions stipulated in Ref. []Sl| , 
but this is not clear to us at this stage. 

As another example, J N can be used as a figure of merit 
in designing optimized quantum control and/or quantum 
error correction strategies: One is typically interested in 
determining a quantum operation 6 that minimizes the 
averaged distance between the elements of a set of noisy 
quantum states p i and a predefined set of target quantum 
states <7i. In this context, it would be interesting to 
investigate if distance measures based on J N would lead 
to any advantage in terms of computation time. Clearly, 
this has potential applications to the implementation of 
real-time quantum technologies. 

Yet another possible direction of research consists of 



employing 3> as a distance measure between quantum 
operations — as opposed to quantum states — via the 
isomorphism between quantum states and CPTP maps 
|p2[ |63| . In this regard, it is worth investigating whether 
distance measures based on 3 r N would satisfy the six cri- 
teria proposed in Ref. . Remarkably, from the results 
of the present work and Ref. [Q, a few strengths of 
iFjv-based measures can already be anticipated. Of spe- 
cial significance are the fulfillment of the criteria "easy 
to calculate" and "easy to measure." Along these lines, 
some operational meaning for $ N would also be highly de- 
sirable. Although we do not presently have a compelling 
physical interpretation of 3^ , it is not inconceivable that 
one can be found in an analogous way to 5" ElJ . 
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APPENDIX A: METRICS 

From a mathematically rigorous viewpoint, a distance 
measure ID on a set S is a function D : S x S — > K such 
that for every a, 6, c £ S the following properties hold. 

(Ml) D(a,b) > (non- negativity), 

(M2) ID (a, b) = iff a = b (identity of indiscernible), 

(M3) D(a,b) = D(b,a) (symmetry), 

(M4) ID (a, c) < ID (a, b) + ID (6, c) (triangle inequality). 

Any such function is called a metric. 

APPENDIX B: PROOFS 



1. Proof of Proposition III. 2 



In this appendix the joint concavity of J N is established 
via the proof of Proposition III. 2. 
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Proof. Differentiating Eq. ( p0| ) twice with respect to x, 
we obtain 

d 2 F(x) 



dx 2 



= 2u ■ v 
+ 2 



g{x) + f(x)- 



dx 2 
df{x) dg(x) 



dx 2 



dx dx ' 

where, for convenience, we define the functions f(x) := 



\Jl — \\r + xu\\ 2 and g(x) := sJT 
some computation we find that 



\s + xv\ 



d 2 F(x) 
dx 2 



where 



3l(x) := 2u ■ v — 



g{x)u 2 f{x)v 2 



f{x) g(x) 

u ■ (f + xu)v ■ (s + xv) 



f(x)g(x) 

g(x)[u ■ (f + xu)] 2 f (x)[v ■ (s + xv)Y 



lf(xW 



[g(xW 



(Bl) 

r) := 
After 

(B2) 
(B3) 

(B4) 



The negative semidefiniteness of d 2 F{x)/dx 2 in the range 
x £ [0, 1] can be observed if $i{x) and #2 0*0 are written 
in the following alternative form: 

2 

ffi(aO 




f(x)g{x 



lf(x] 



«■(?' + OT TT V ' ( S + 

9(x) 



□ 



The supermultiplicative property is obtained by show- 
ing the positive semidefiniteness of the right-hand side of 
the above expression. This is the content of the following 
proposition. 

Proposition B.l. For < a,b,c,d < 1, we have 



y/(l - ab)(l - cd) > y/(l - a)(l - b)(l - c)(l - d) 

+ v/ac(l - b)(l - d) + y/bd(l-a)(l-c). (B5) 

Proof. First note that if any of the variables equals 1, 
then the validity of the inequality i s im mediate. For 
example, let d = 1 so that inequality (B5) reduces to 



y/(l-ab){l-c) > y/b(l-a)(l-c). 



(B6) 



This is trivially satisfied for all < a, b, c < 1. In what 
follows, we restrict ourselves to < a,b,c,d < 1 and 
show that inequality (|B5| ) is equivalent to the standard 
inequality of arithmetic and geometric means (hereafter 
referred as the AM-GM inequality). This inequality is 
just an expression of the fact that the geometric mean of 
a list of non-negative real numbers is never larger than 
the corresponding arithmetic mean. 

Apply the substitution a' = 1 — a (similarly for b' , d, 
and d'; note that < a',b',c',d' < 1) to inequality (Jb|) 
and divide the result by v a'b'c'd' to get the equivalent 
inequality 



y / {l + A + B)(l + C + D) > 1 + VIC + VBD, (B7) 

where we have defined A = 1/af — 1 (similarly for B, C, 
and D; note that < A,B,C,D < oo). Squaring the 
inequality above we find 

A+C B + D AD + BC 



> VAC- 



2. Proof of supermultiplicativity of 5F A 



2 2 2 

which is clearly a sum of three AM-GM inequalities. □ 



r BD+VABCD 
(B8) 



To prove thctt U'jv 

is supermultiplicative, we first define 
Ti ■= Tr(pf) and s, := Tr(crf), such that < r it Si < 1 
(note that here we use r, instead of rf as the norm square 
of r*i, likewise for Sj). Straightforward algebra gives 

3n(Pi <8> P 2 , cr i ® &i) - 3 r N (p 1 ,cr 1 )3 : N (p 27 cr 2 ) = 



VAl - nr 2 )(l - sis 2 ) 

- v / (l-ri)(l-s 1 )(l-r 2 )(l-s 2 ) 

-Tv(p 1 cr 1 )^(l-r 2 )(l- S2 ) 

- Tr(p 2 o- 2 )y(i -n)(l- Sl ). 

A direct application of Cauchy-Schwarz's inequality 
Tr^crj) < ^/rjS, gives 

9>(Pl ® P 2 ,CTl ® 0-2) - J N {p 1 ,(Tl)'J N (p 2 ,(T2) > 

- v /(l-r 1 )(l- Sl )(l-r 2 )(l-s 2 ) 



V r l s l(l - r 2)(l - « 2 ) - V r 2S2(l - 1"i)(l - Si). 



3. Proof of the metric property of B[J] and C[3"] 

In the following, we give an alternative demonstration 
of the metric properties of B[3] and C[5] (see Refs. |3^, 
Q for the standard proofs). Our proof consists of a 
simple application of Theorem [III.1 due to Schoenberg. 

Proposition B.2. The functionals B[3 r (p, cr)] and 
C[2(p,tr)], defined in Eq. (|) and Eq. (g), are metrics 
for the space of density matrices. 

Proof. Let K[3 r (p, cr)] represent either B[3(p,cr)} or 
C[$(p, cr)] for brevity. As with 5F, it is easy to check 
that K 2 [^(p, cr)] is symmetric in its two arguments and 
that K 2 [3(p, cr)] > with saturation iff p = cr. So, ac- 



cording to Theorem III.l| , K[3 r (p, cr)] is a metric if for any 
set of density matrices {p.j}f =1 (n > 2) and real numbers 
{ c i}"=i such that Yl7=i Ci ~ ^' ^ IS ^ mc ^ na ^ 



K 2 [7{p llPj )] Cl c 3 



< 0. 



(B9) 
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To prove this, we derive an upper bound for 
K 2 [3 r (p i , pj)], which can be easily seen to satisfy the con- 
dition above. First, note that 



Tr 



(VpI 



= Tr 



> 



Tr 



/Piy/Pj y 



= A(p t , Pj ), (BIO) 



where the first equality follows from the definition |A| := 
\/AtA for every matrix A and the inequality from the 
fact that Tr(|A|) = maxu | Tr(UA)| (the maximization 
runs over unitary matrices U f34|). Then, it follows 
that 

B 2 [J] = 2(1 - V?) < 2(1 - 3) < 2(1 - A), (Bll) 
C 2 [9] = l-3 r <l-A<2{l-A), (B12) 

or, in our more compact notation, A' 2 [5F] < 2(1 — A). 

Now, replacing K 2 [3 r (p i , p -)] with the above upper 
bound in the left-hand side of Eq. (|l39| ), it is easy to 
obtain the desired inequality: 



£,.7=1 



2 - 2Tr 



)]}., 



2Tr 



E 1 



< 0, (B13) 



where the equality is obtained by using the fact that 
ST=i c i ~ 0, the linearity of the trace operation, and 
the hermiticity of Ci-^fpl® J~f>i- □ 

Finally, let us just mention that besides establishing 
the metric properties of -B[5] and C[3] : the present 



proof also establishes y 2 — 2 [Tr (-y/p-y/ir)] as a me- 
tric for the space of density matrices. In fact, by a si- 
milar application of Schoenberg's theorem, the quantity 

H(p 1 cr) := y^2 — 2 Tr (yfpy/a^) can also be shown to be 
a metric. 



APPENDIX C: MATLAB CODES 

In this appendix, we present the Matlab codes that we 
have used to compute the various functions involved in 



the numerical experiment presented in Sec. [V 
For rho and sigma density matrices, 

• 3> was computed using 

Fn = real( rho ( : ) ' *sigma( : ) ... 

+ sqrt((l - rho(:) '*rho(:))* 
(1 - sigma ( : ) '*sigma( : ) )) ); 

• 5" was computed using 



[V, D] = eig(rho); 

sqrtRho = V*diag(sqrt(diag(D)))*V ; 
F = sum( sqrt(eig(Hermitize( ... 

sqrtRho*sigma*sqrtRho) ) ) )~2; 

Here sqrtRho*sigma*sqrtRho is not quite Hermi- 
tian due to small numerical errors. We therefore 
employ the function Hermitize(M) = (M+M' )/2 to 
turn the almost-Hcrmitian matrix into a Hermitian 
one this causes Matlab to select a more efficient 
algorithm for the diagonalization. 

• D was computed using 

D=0 . 5*sum(abs ( eig(rho-sigma) )); 

• Q was computed using 

[Vr ,Drho] =eig(rho) ; Dr=diag(Drho) ; 

[Vs ,Dsigma] =eig(sigma) ; Ds=diag(Dsigma) ; 

A = abs(Vr'*Vs) . "2; 

[x,Q]=fminbnd(@(s) ... 

(Dr. ' .~s)*A*(Ds.~(l-s)), 0, 1); 

The algorithm used here follows from the formula 
for Tr(p s cr 1 ~ s ) given in the section entitled conve- 
xity in s of Ref . |GM . 
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